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Abstract 

The Hong-Ou-Mandel (HOM) dip plays an important role in recent linear optics 
experiments. It is crucial for quantum computing with photons and can be used to 
characterise the quality of single photon sources and linear optics setups. In this 
paper, we consider generalised HOM experiments with N bosons or fermions passing 
simultaneously through a symmetric Bell multiport beam splitter. It is shown that for 
even numbers of bosons, the HOM dip occurs naturally in the coincidence detection 
in the output ports. In contrast, fermions always leave the setup separately exhibiting 
perfect coincidence detection. Our results can be used to verify or employ the quantum 
statistics of particles experimentally. 

1 Introduction 

The two-photon Hong-Ou-Mandel (HOM) dip has been demonstrated first in 1987 |1J. In 
their experiment, Hong, Ou and Mandel sent two identical photons simultaneously through 
the separate input ports of a 50 : 50 beam splitter. Each output port contained a photon 
detector. Surprisingly, no coincidence detections within the temporal coherence length of 
the photons, i.e. no simultaneous clicks in both detectors, were recorded. Crucial for the 
observation of this effect was the identicalness of the pure quantum states of the input 
photons, which differed only in the directions of their wave vectors. This allowed the photons 
to become indistinguishable and to interfere within the setup. The detectors could not resolve 
the origin of each observed photon. 

Due to the nature of this experiment, the HOM dip was soon employed for quantum 
mechanical tests of local realism |2} and for the generation of postselected entanglement 
between two photons jH] • Linear optics Bell measurements on photon pairs rely intrinsically 
on the HOM dip jUEUHl, which has also been a building block for the implementation of linear 
optics gates for quantum information processing with photonic qubits [7|. Shor's factorisation 



1 



N 



Figure 1: (a) Pyramidal construction of an N x N Bell multiport beam splitters for photons 
using a network of beam splitters and phase plates, (b) Alternatively, an N x N multiport 
can be constructed out of iV fibres spliced together. 

algorithm jS] , for example, relies on multiple path interference to achieve massive parallelism 
jH] and multiphoton interference has to play a crucial role in any implementation of this 
algorithm using linear optics. 

Since it requires temporal and spatial mode-matched photons, observing the HOM dip 
for two photons is also a good test of their indistinguishability. HOM interference has been 
applied to characterise recently introduced sources for the generation of single photons on 
demand by testing the identicalness of successively generated photons [101 El Ej . Another 
interesting test based on the HOM dip has been studied by Bose and Home, who showed 
that it can reveal whether the statistics of two identical particles passing through a 50 : 50 
beam splitter is fermionic or bosonic [T3] . 

Motivated by the variety of possible applications of the two-photon HOM dip, this paper 
investigates generalised HOM experiments. We consider a straightforward generalisation of 
the scattering of two particles through a 50 : 50 beam splitter, namely the scattering of 
iV particles through a symmetric N x N Bell multiport beam splitter. While numerous 
studies on N photon interference in the constructive sense, i.e. resulting in the enhancement 
of a certain photon detection syndrome, have been made (see e.g. Refs. 0), not much 
attention has been paid to multiple path interference in the destructive sense. Recently 
Walborn et al. studied so-called multimode HOM effects for photon pairs with several inner 
degrees of freedom, including the spatial and the polarisation degrees of freedom [TH Ho"]. 
A notable example for destructive HOM interference has been given by Campos |16j . who 
studied certain triple coincidences in the output ports of an asymmetric 3x3 multiport 
beam splitter, which is also known SIS db fritter. 

We consider bosons as well as the simultaneous scattering of fermions. The difference 
between both classes of particles is most elegantly summarised in the following commutation 
rules. While the annihilation and creation operators cij and a\ for a boson in mode i obey 
the relation 

[cii, Oj] = diOj — alai = 5ij and [a\, al] = [a*, a,-] =0 V i, j (1) 
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with 5ij = for i ^ j and 5u = 1, the annihilation and creation operators a, and a\ of 
fermionic particles obey the anticommutation relation 

{oj, Oj} = CLiOj + a]oj = <5y and {af , aj} = {aj, a,} = V i, j . (2) 

Here i and j refer to the inner degrees of freedom of the particles, like their respective path, 
polarisation, spin, frequency or energy. 

Possible realisations of the N x N Bell multiport for photons, behaving like bosons, are 
shown in Figure They may consist of a network of beam splitters and phase plates [T71ITH] 
but can also be made by splicing N optical fibers [THJ. Such spliced fibre constructions are 
commercially available and can include between 3 and 30 input and output ports. The main 
feature of the symmetric Bell multiport is that a photon entering any of the input ports is 
equally likely redirected to any of the possible output ports. It can therefore be used to 
produce higher dimensional EPR correlations JB] and for the generation of so-called NOON 
states with special applications in lithography 1201 El 122] • Moreover, the N x N Bell 
multiport can be used to prepare a great variety of multiphoton entangled states (221 1211 
and to maximise the success probability when teleporting photonic qubits jjj- 

Furthermore, Bell multiport beam splitters exist for a wide variety of fermions and bosons. 
For example, multiports for bosonic or fermionic atoms can consist of a network of electrode 
wave guide beam splitters on an atom chip [23]. Multiports for electrons, who behave like 
fermions, can be realised by fabricating a network of quantum point contacts |2S], split gates 
|2*7j or semiconductor multiterminal nanostructures (2H1 1211 acting as two-electron beam 
splitters and can be used, for example, to detect two-electron Bell states [3U|. Specially 
doped optical fibres have recently been introduced in the literature and are expected to 
constitute beam splitters for "fermion-like" photons [31] . 

As in the original HOM experiment [Tj, we assume in the following, that a particle 
detector is placed in each output port of the scattering beam splitter array. The incoming 
particles should enter the different input ports more or less simultaneously and such that 
there is one particle per input port. Moreover, we assume that the particles are identical. It 
is shown, that it is impossible to observe a particle in each output port for even numbers N 
of bosons. We denote this effect of zero coincidence detection as the generalised HOM dip. It 
is also shown that fermions always leave the setup separately exhibiting perfect coincidence 
detection. Since the interference behaviour of both types of particles is very different, the 
Bell multiport can be used to reveal their quantum statistics. 

This paper is organised as follows. In Section |21 we introduce the theoretical description 
of particle scattering through a symmetric Bell multiport. Section El describes the scattering 
of two particles through a 50 : 50 beam splitter as an example. In Section^ we derive the 
condition for the generalised HOM dip for bosons and analyse the scattering of fermions 
through the same setup for comparison. Finally we conclude our results in Sectional 
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2 Scattering through a Bell multiport beam splitter 



Let us first introduce the notation for the description of the scattering of N particles through 
a passive setup. We proceed in close analogy to the scattering of photons through a linear 
optics network [T5J Suppose each input port i is entered by a particle with creation 
operator a\. Then the input state of the system equals 



N 



) = IHlO>. ( 3 ) 



where |0) is the vacuum state with no particles in the setup. Moreover, S denotes the unitary 
scattering matrix, which connects the input state to its output state 

|<M = S\<f) in }. (4) 

Using Eq. Q and the relation SrS = 1 then yields 

N 

|0 out ) = Sal&Sak- ...■&Sal f S*S\0) = ]l Sat& \0), (5) 

i=i 

since S\0) = |0) for the considered passive setup. As long as no particles enter the system, 
it remains in its vacuum state. 

The Bell multiport beam splitter directs each incoming particle with equal probability 
to all output ports. To describe this, we introduce the matrix elements Uji, which give the 
amplitude for the distribution of a particle from input port i to output port j. Especially, 
for anJVx JV Bell multiport one has 



where ujn denotes the iV-th root of unity 

uj n = e 2i ^ N . (7) 

The corresponding matrix U is also known as the discrete Fourier transform matrix and has 
been widely considered in the literature [THl UHl 123 E2 ■ Proceeding as in Section II. D of 
Ref. [TSj . it can easily be verified that U is unitary as well as symmetric. 

If bj denotes the creation operator for a single particle in output port j, the definition of 
the scattering matrix U implies that 

Sa\S^ = Y.U ]l b). (8) 

3 

Inserting this into Eq. ©, we obtain 

N I N \ 

i^out) = n n^b) no). (9) 

1=1 \ ,=1 / 
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Figure 2: (a) HOM dip for two bosons scattering through a 50 : 50 beam splitter, (b) Perfect 
coincidence in the output ports for fermion scattering. 

This equation describes the independent redirection of the incoming particles to their poten- 
tial output ports. Conservation of the norm of the state vector is provided by the unitarity of 
the transition matrix U. Note that up to now, we have not invoked any assumptions about 
the nature of the particles. The formalism in this section applies to bosons and fermions 
equally. 



3 HOM interference of two particles 

Before analysing the general case, we motivate our discussion by considering two identical 
particles entering the different input ports of a 50 : 50 beam splitter. For N = 2, the 
transition matrix (JHJ) becomes the Hadamard matrix 



U = ±A (io) 




and the input state © becomes \(f)- m ) = a\a\ |0). Note that local measurements on this 
input state cannot reveal any information about the bosonic or fermionic nature of the two 
particles. However, using Eq. (JUJ), we find that the beam splitter prepares the system in the 
state 

i^out) = \ ip\ + bt)(b\ - bt) io> = \ [(b\) 2 - b\b\ + b\b\ - (b\f\ io> . (ii) 

This state no longer contains any information about the origin of the particles, since each 
incoming one is equally likely transferred to any of the two output ports. Passing through 
the setup, the input particles become indistinguishable by detection (see Figure EJ)- Their 
quantum statistics can now be revealed using local measurements. 

Bosons obey the commutation law (Q). Using this, the output state becomes 



'■>OVLt 



) = \ W - (4) 2 1 |0> , (12) 



which implies a zero-coincidence count rate at the output ports. The particles bunch together 
in the same output port and exhibit the famous HOM dip (see Figure |2fa)). In contrast, 
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(13) 



implies perfect particle coincidence. This means that the fermions always arrive in separate 
output ports and never bunch together (see Figure 0(b)) . A 50 : 50 beam splitter can 
therefore indeed be used to distinguish bosons and fermions [T3] . 



We now consider the general case of N particles passing through aniVx JV Bell multiport 
beam splitter. As in the N = 2 case, the setup redirects each incoming particle with equal 
probability to any of the possible output ports, thereby erasing the information about the 
origin of each particle and making them indistinguishable by detection. For even numbers of 
bosons, this results in the generalised HOM dip and zero coincidence detection. In contrast, 
fermions leave the setup always separately, thus demonstrating maximum coincidence de- 
tection. Observing this extreme behaviour can be used, for example, to verify the quantum 
statistics of many particles experimentally. 

4.1 Bosonic particles 

In order to derive the necessary condition for the appearance of the generalised HOM dip for 
even numbers of bosons, we calculate the output state of the system under the condition 
of the collection of one particle per output port. Each term contributing to the projected 
conditional output state \4> pro ) can be characterised by a certain permutation, which maps 
the particles in the input ports 1, 2, iV to the output ports 1, 2, N. In the following, 
we denote any of the AH permutations by a with a(i) being the z-th element of the list 
obtained when applying the permutation a onto the list {1, 2, iV}. Using this notation, 
|0 P ro) equals up to normalisation 



is the probability to detect one particle per output port. It is therefore also the probability 
for observing coincidence counts in all iV detectors. 

Up to now, the nature of the particles has not yet been taken into account. Using the 
commutation relation ([!]) for bosons, the conditional output state (fH| becomes 



4 Multiparticle HOM interference 




(14) 



The norm of this state has been chosen such that 




(15) 



N 




(16) 



i=i 
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with the permanent of the matrix U defined as 

N 

permC/ = ^n^«- ( 17 ) 

V 1 = 1 

The permanent of a matrix is superficially similar to the determinant. However, there exist 
hardly any mathematical theorems that can simplify the calculation of the permanent of an 
arbitrary matrix. 

To derive a condition for the impossibility of coincidence detections, we have to see when 
the probability (|T5j) equals zero. Using Eq. (|T5|). we find 

-Pcoinc = I perm £7 | 2 . (18) 

The key to the following proof is to show that the transition matrix U of the Bell multiport 
possesses a certain symmetry such that its permanent vanishes in certain cases. Suppose the 
matrix U is multiplied by a diagonal matrix A with matrix elements 

A jk = w j x 8 jk . (19) 

This generates a matrix AU with 

N 

(AC/),, = £ A jk U ki = A n U 3l = jft ut l)t ■ (20) 

k=l 

We now introduce the modulus function defined as mod^x) = j, if x — j is dividable by N 
and < j < N. Since = u% = 1, the matrix elements (|2Uj) can be expressed as 

(\TT\ — 1 , ,0'-l)(mod JV (i)+l-l) / -,s 
{AU)ji — -j=U N . [Zi) 

Note that the function a(i) = modTv(^) + 1 maps each element of the list {1, 2, ...N — 1, iV} 
respectively to the list {2, 3, ...N, 1}. A comparison with Eq. JHJ) therefore shows that 

{AU)ji — Uja(i) ■ (22) 

In other words, the multiplication with A amounts to nothing more than a cyclic permutation 
of the columns of the matrix U. Taking the cyclic permutation symmetry of the permanent 
of a matrix (see definition ()17|) ) into account, we obtain 

perm U = perm (AU) . (23) 

However, we also have the relation 

perm (AU) = perm A ■ perm U (24) 
with the permanent of the diagonal matrix A given by 



perm A = g & = = < (K+ " /2 = e"<~«> = { J | ™$ ^ (25) 



For TV being even, a comparison of Eqs. (|23[) - (|25|) reveals that 

perm = —perm U = . (26) 

As a consequence, Eq. fITHj) implies that P C omc = 0. Coincidence detection in all output 
ports of the setup is impossible for even numbers of bosons. This is not necessarily so, if the 
number of particles is odd. For example, for N = 3 one can check that there is no HOM 
dip by calculating perm U explicitly. Campos showed that observing a HOM dip for N = 3 
is nevertheless possible with the help of a specially designed asymmetric multiport beam 
splitter [To] . 

4.2 Fermionic particles 

Fermions scattering through a Bell multiport beam splitter show another extreme behaviour. 
Independent of the number N of particles involved, they always leave the setup via different 
output ports, thereby guaranteeing perfect coincidence detection. As expected, particles 
obeying the quantum statistics of fermions cannot populate the same mode. 

Again, we assume that each input port is simultaneously entered by one particle and 
denote the creation operator of a fermion in output port i by b\. Proceeding as in Section 
14.11 one finds again that the output state of the system under the condition of the collection 
of one particle per output port is given by Eq. (JHJ). To simplify this equation, we now 
introduce the sign function of a permutation with sgn(<r) = ±1, depending on whether the 
permutation a is even or odd. An even (odd) permutation is one, that can be decomposed 
into an even (odd) number of integer interchanges. Using this notation and taking the 
anticommutator relation for fermions (J2J) into account, we find 

|0 pro ) = £sgm»( l[U a{i)i bi) |0> . (27) 

cr \i=l ) 

A closer look at this equation shows that the amplitude of this state relates to the determinant 
of the transformation matrix given by 

TV 

dettf = $>gn(cr)n^(i)i- (28) 

cr i=l 

Since U is unitary, one has |det U\ — 1 and therefore also, as Eq. (|T5)l shows, 

Pcoinc = I det £/T = 1 . (29) 

This means, that fermions always leave the system separately, i.e. with one particle per 
output port. In the above, we only used the unitarity of the transition matrix U but not 
it's concrete form. Perfect coincidence detection therefore applies to any situation where 
fermions pass through an N x N multiport, i.e. independent of its realisation. 
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5 Conclusions 



We analysed a situation, where N particles enter the N different input ports of a symmetric 
Bell multiport beam splitter simultaneously. If these particles obey a fermionic quantum 
statistics, they always leave the setup independently with one particle per output port. This 
results in perfect coincidence detection, if detectors are placed in the output ports of the 
setup. In contrast to this, even numbers N of bosons have been shown to never leave the 
setup with one particle per output port. This constitutes a generalisation of the two-photon 
HOM dip to the case of arbitrary even numbers N of bosons. The generalised HOM dip is 
in general not observable when iV is odd. 

The proof exploits the cyclic symmetry of the setup. We related the coincidence detection 
in the output ports to the permanent or the determinant of the transition matrix U describ- 
ing the multiport, depending on the bosonic or fermionic nature of the scattered particles. 
Although the definition of the permanent of a matrix resembles that of the determinant, 
there exist almost no theorems to simplify their calculation I35j. In fact, the computa- 
tion of the permanent is an NP-complete problem and much harder than the calculation of 
a determinant, which is only a P problem in complexity. Experimental setups involving the 
scattering of bosons through a multiport therefore have important applications in quantum 
information processing. 

For example, part of the linear optics quantum computing scheme by Knill, Laflamme 
and Milburn ;7J is based on photon scattering through a Bell multiport beam splitter. In 
contrast to this, the scattering of non-interacting fermions through the same corresponding 
circuit, can be efficiently simulated on a classical computer P^l I37j. Moreover, the quan- 
tum statistics of particles has been used for a variety of quantum information processing 
tasks such as entanglement concentration [3H] and entanglement transfer Completely 
new perspectives might open when using setups that can change the quantum behaviour of 
particles and convert, for example, photons into fermions [23] . 

Finally, we remark that observing HOM interference of many particles is experimentally 
very robust. Our results can therefore also be used to verify the quantum statistics of 
particles experimentally as well as to characterise or align an experimental setup. Testing the 
predicted results does not require phase stability in the input or output ports nor detectors 
with maximum efficiency The reason is that any phase factor that a particle accumulates in 
any of the input or output ports contributes at most to an overall phase factor of the output 
state |0out)- However, the coincidence statistics is sensitive to the phase factors accumulated 
inside the multiport beam splitter as they affect the form of the transition matrix (jUJ). The 
case, where the input ports of the setup are not entered by perfect one-particle states but 
by mixtures containing also a vacuum component can be analysed, in principle, using the 
methods introduced by Berry et al. in Refs. j4(J[ I41j. 
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